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ABSTRACT 

We model the spectra and surface brightness distributions for the Lya radiation expected from protogalaxies 
that are caught in the early stages of their assembly. We use the results of a companion paper to characterize 
the radiation emerging from spherically collapsing gas clouds. We then modify these spectra to incorporate the 
effect of subsequent resonant scattering in the intergalactic medium (IGM). Using these models, we interpret a 
number of recent observations of extended Lya blobs (LABs) at high redshift. We suggest, based on the angular 
size, energetics, as well as the relatively shallow surface brightness profiles, and double-peaked spectra, that 
several of these LABs may be associated with collapsing protogalaxies. We suggest two follow-up observations to 
diagnose the presence of gas infall. High S/N spectra of LABs should reveal a preferential flattening of the surface 
brightness profile at the red side of the line. Complementary imaging of the blobs at redshifted Ha wavelengths 
should reveal the intrinsic Lya emissivity and allow its separation from radiative transfer effects. We show that 
Lya scattering by infalling gas can reproduce the observed spectrum of Steidel et al's LAB2 as accurately as a 
recently proposed outflow model. Finally, we find similar evidence for infall in the spectra of point-like Lya 
emitters. The presence of scattering by the infalling gas implies that the intrinsic Lya luminosities, and derived 
quantities, such as the star-formation rate, in these objects may have been underestimated by about an order of 
magnitude. 

Subject headings: galaxies: formation - galaxies: halos - quasars: general - radiative transfer - cosmology: 
theory - intergalactic medium 



1. INTRODUCTION 

Copious Lya emission may be associated with the early stages 
of the formation of galaxies. The epoch before the onset of sig- 
nificant star formation in the history of an individual galaxy, 
hereafter referred to as the 'protogalactic' stage, may be ac- 
companied by the release of the gravitational binding energy 
through Lya cooling radiation (Haiman et al. 2000; Fardal et al. 
2001). The extended Lya "fuzz" surrounding these protogalax- 
ies may be strongly enhanced in the presence of a central ion- 
izing source, such as a quasar (Haiman & Rees 2001; Haiman 
2004), as evidenced by a handful of recent observations (Bunker 
et al. 2003a; Weidinger et al. 2004, 2005; Christensen et al. 
2006). As the number of known high redshift Lya emitters, 
both point-like (e.g. Rhoads et al. 2000; Rhoads & Malhotra 
2001; Ajiki et al. 2003; Kodaira 2003; Hu et al. 2004; Rhoads 
et al. 2004; Ouchi 2005) and extended (e.g Steidel et al. 2000; 
Matsuda 2004; Wilman et al. 2005; Dey et al. 2005; Nilsson 
et al. 2005; Matsuda et al. 2006) grows, there is an increasing 
need for diagnosing the presence of any large-scale gas infall. 

In Dijkstra et al. (2006, hereafter paper I) we presented Monte 
Carlo simulations of the transfer of Lya photons through spher- 
ically symmetric models of neutral collapsing gas clouds. These 
models are intended as simplified descriptions of protogalax- 
ies in the process of their assembly. We computed the spec- 
tra and surface brightness profiles emerging from clouds with 
a range of properties. For clouds without any associated ion- 
izing source, generating Lya photons through cooling radia- 
tion alone, we explored two extreme models. In the 'extended 
models', the gas continuously emits a significant fraction of 
its gravitational binding energy in the Lya line as it navigates 
down the dark matter potential well. In the "central" models, 
Lya photons are generated only at the center of the collapsing 
cloud. We then studied how the results of these calculations 



were affected by the presence of a central ionizing source. The 
main results that are relevant to the present paper can be sum- 
marized as follows: 

(i) In the absence of an ionizing source (i.e. stars or a quasar) 
the emergent Lya spectrum is typically double-peaked and asym- 
metric. Transfer of energy from the collapsing gas to the Lya 
photons -together with a reduced escape probability for photons 
in the red wing- causes the blue peak to be significantly en- 
hanced, which results in an effective blueshift of the Lya line. 
This blueshift may easily be as large as ~ 2000 km s" 1 . The 
total blueshift of the line increases with the total column den- 
sity of neutral hydrogen of the gas and with its infall speed. The 
prominence of the red peak is enhanced towards higher redshift, 
lower mass and lower infall speeds. If detected, the shift of the 
red peak relative to the line center potentially measures the gas 
infall speed. 

(ii) In the "central" models, the Lya radiation emerges with 
a steeper surface brightness profile, and with a larger blueshift, 
than in the "extended" models. 

(iii) A steepening of the surface brightness distribution to- 
wards bluer wavelengths within the Lya line is indicative of 
gas infall. 

(iv) In clouds that are fully ionized by an embedded ioniz- 
ing source, the effective Lya optical depth is reduced typically 
to < 10 3 . In these cases, the Lya emission emerges with a 
nearly symmetric profile, and its FWHM is determined solely 
by the bulk velocity field of the gas. Its overall blueshift is re- 
duced significantly (and may vanish completely for very bright 
ionizing sources). As a result, infall and outflow models can 
produce nearly identical spectra and surface brightness distri- 
butions, and are difficult to distinguish from one another. 

For a detailed description of several other features of the 
emergent Lya radiation the reader is referred to paper I. The 
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models discussed in paper I are oversimplified, for a discussion 
on the uncertainties this may introduce the reader is referred to 
paper I, specifically section 8.2. The results presented in paper I 
generally ignored the transfer through the intergalactic medium 
(IGM). Since the IGM is optically thick at the Lya frequency 
at z 3, to allow comparison of these results with the observa- 
tions, the scattering of Lya photons through the IGM needs to 
be modeled. Note that we define the IGM to commence beyond 
the virial radius of an object. The intrinsic spectrum refers to 
the spectrum as it emerges at the virial radius. In the present 
paper, we adopt a simple model of the IGM gas density and 
velocity distribution expected near protogalaxies, based on the 
models of Barkana (2004). Combining the results obtained in 
paper I with the scattering in the IGM, we compare our mod- 
els to several observations in the literature. We focus on inter- 
preting those extended Lya blobs (LABs) that have no contin- 
uum source, or surround a radio quiet quasar. While LABs are 
known to exist around radio loud quasars, they are attributable 
to known jets and supernova-driven outflows, and we do not 
attempt to model them here (e.g. Hu et al. 1991; Wilman et al. 
2000). We also compare our results with the recently discov- 
ered high redshift (z > 5) Lya emitters, which generally appear 
point-like. Despite the over-simplified nature of our model of 
the protogalactic gas and the IGM, we successfully reproduce 
several observed Lya spectral features. 

The rest of this paper is organized as follows. In § 2, we 
describe our model of the IGM and discuss its impact on our 
results from paper I. In § 3, we compare our model predictions 
with the properties of several known high-redshift Lya emit- 
ting sources. In § 4, we discuss our results and their cosmo- 
logical relevance. We also discuss the possibility of imaging 
the protogalaxies in the Ha. Since this line should not suffer 
from attenuation in the IGM, it may have a surface brightness 
comparable to (or even exceeding) that of the Lya line. Fi- 
nally, in § 5, we present our conclusions and summarize the 
implications of this work. The parameters for the background 
cosmology used throughout this paper are Q m = 0.3, = 0.7, 
Ct b = 0.044, h = 0.7, based on Spergel et al. (2003). 

2. THE IMPACT OF THE IGM 

2.1. Impact of a Static IGM 

In the simplest model, employed occasionally in paper I, Lya 
photons are injected into a statistically uniform IGM, following 
Hubble expansion. Photons with an original wavelength blue- 
ward of the Lya line center redshift into resonance and may get 
scattered out of our line of sight. Lya photons that are initially 
redward of the line center are not scattered resonantly (scatter- 
ing in the red damping wing is not important in the present con- 
text, at redshifts at which the IGM is highly ionized). We fol- 
low paper I and define x as a dimensionless frequency variable, 
which measures the offset from the line center in units of the 
Doppler width, x = (v-vq)/ l\v D . Note that positive/negative 
values of x correspond to blueshifted/redshifted Lya photons. 
The IGM then allows a transmission of a fraction (exp (— r)) and 
1 with x > and x < of photons emerging from the collaps- 
ing cloud, respectively. Here, (exp(--r)} denotes the mean IGM 
transmission as it has been determined observationally in stud- 
ies of the Lya forest. Table 1 shows (exp(-r)) as a function of 
redshift. The second row shows ay, which is the reported mea- 
surement error (1 —a) in /. This Table shows that attenuation 
by the IGM is significant at z k 3. 



Table 1 
Transmission of IGM 

_z 2.4 1 3.0' 3.8' 4.5 2 5.2 2 5.7 2 6.0 2 

/ = (exp(-r)) 08 068 049 021 O09 O07 O004 

oy 0.02 0.003 0.003 

(1) From McDonald et al. (2001) 

(2) From Songaila & Cowie (2002); Fan et al. (2002) 
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FIG. 1. — This drawing illustrates the effect of the impact parameter on 
the IGM transmission. The IGM is assumed to undergo spherically symmet- 
ric infall. In this case, photons emerging from the edge of our object (in the 
direction along the line of sight) see a projected version of the IGM velocity 
field. 



2.2. Impact of a Dynamic IGM 

Protogalaxies are likely to reside in relatively massive halos, 
corresponding to high-er density peaks in the IGM. The gravita- 
tional field of the halos causes material on large scales to decou- 
ple from the Hubble flow and fall in toward the halo. Barkana 
(2004) has modelled the density and velocity distribution of the 
IGM around nonlinear halos. In these models, around the mas- 
sive halos typically considered in paper I, infall may occur up 
to about ~ 5 virial radii. Photons that escape the protogalaxy 
with a small redshift, x <^ 0, may then be resonantly scattered by 
the infalling gas. As a result, significant attenuation can extend 
into the red part of the observed spectrum. To give a quantita- 
tive example (see Figure 2 and 4 in Barkana 2004), for a dark 
matter halo with a total mass of 10 12 M Q , the overdensity of the 
IGM is expected to be S ~ 20- 30 at the virial radius, and still 
5^2-3 at ~ 10 virial radii from the object. The infall velocity 
is equal to the circular velocity at the virial radius. The IGM 
is at rest with respect to the virialized object approximately 5 
virial radii away. The radius at which the peculiar velocity of 
the IGM vanishes (and is therefore comoving with the Hubble 
flow ) is further out, and is mass dependent. 

The impact of the IGM will be further modified by the prox- 
imity effect near bright QSOs. For bright QSOs, the size of 
the local HII region can extend beyond <~ 5R v \ r , significantly 
reducing resonant absorption from the dynamically biased part 
of the IGM (Cen & Haiman 2000; Madau & Rees 2000). In our 
models below, we allow for the presence of such a cosmolog- 
ical HII region. Finally, the impact of the IGM can also vary 
from location to location, due to stochastic density and ionizing 
background fluctuations. We will restrict our study to study the 
mean effect of the IGM, and not investigate such variations in 
the present paper. 

2.3. Modelling the Dynamic IGM 

To illustrate the impact of the IGM, we model the transmis- 
sion around a halo with a total mass M tot = 5.2 x 10 11 M (cor- 
responding to the mass of fiducial model I. in paper I), in the 
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presence or absence of a bright quasar. Note that for increasing 
masses, the IGM would have increasing peculiar velocities, and 
would affect a wider range of frequencies, extending further to 
the red side of the line. We assume that the IGM is falling in 
with a velocity profile Vigm(/)> an d that the total number den- 
sity of hydrogen atoms (neutral+ionized) is «igm('")- The opti- 
cal depth for a Lya photon that enters the IGM at frequency x 
is given by 



Red 



z = 3 



z = 5.2 



Blue 



ds higmO) x m (s) cr <f>(x[s]), 



(1) 



where jchi is the neutral fraction of hydrogen by number. We 
take density and velocity profiles which follow the curves shown 
in Figures 2 and 4 in Barkana (2004). We assume a gas temper- 
ature of T = 10 4 ° K (as in paper I, appropriate for photoionized 
gas). For an extended, spherically symmetric source, Lya pho- 
tons propagating along the line of sight "see" a velocity compo- 
nent of the IGM that depends on their impact parameter. This is 
illustrated in Figure 1 . The angle 6> v ; r is the angle on the sky ex- 
tended by the virial radius r v ; r . We adopt a standard expression 
for the virial radius, as a function of halo mass and redshift, 
from Barkana & Loeb (2001). Photons that were received from 
the center of the object, 9 = 0, see the full magnitude of the IGM 
peculiar velocity vector. However, photons that were received 
from \9\ > see only the projection of the infall velocity along 
the line of sight. As a result, the IGM transmission depends on 




t x (9) = (T / dr 



r «igmW Xn(r) 



VigmWt, Ft 1/2 
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Vth 
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where y = 9 x [r v i r /# v ir], and the integral over the line of sight 
s has been converted to an integral over the radial coordinate 
r measuring the distance from the center of the halo to points 
along the line of sight. 

The final step in the calculation is to obtain the neutral frac- 
tion xni(r). As mentioned above, the mean transmission (exp(-r)) 
is known observationally. Here the brackets denote averaging 
over the probability distribution of overdensities, arising from 
small-scale density fluctuations in the IGM. The averaging is 
over a spatial scale corresponding to the spectral resolution. We 
make the simple Ansatz that (exp(-r)) equals exp(-ro), where 
ro is the transmittivity at the mean density and x = 0, 

("lGM)*HI.0VthO'0 



TO : 



H(z) 



(3) 



where in the last step we used ds = v±dx/H(z) (see Appendix 
A of paper I). From this Ansatz, we compute x HI o, the neu- 
tral fraction at the mean density, and, assuming photoioniza- 
tion equilibrium, we compute the background flux J\, g . For gas 
at overdensity 6^0, the neutral fraction then simply scales as 
xui(5) oc (1 + S)J^. In the absence of a quasar, J\, g is constant, 
and in the presence of a quasar, we add the quasar's flux, which 
scales as Jq oc r~ 2 . For a realistic density probability distribu- 
tion function (PDF), the transmission at the mean density can 
differ appreciably from (exp(-r)). However, our simple treat- 
ment, by construction, will give the correct mean transmittiv- 
ity, and effectively assumes only that the gas clumping factor 
C = (p 2 )/{p) 2 due to small-scale inhomogeneities is indepen- 
dent of the overdensity 5 on larger scales (corresponding to a 
spectral resolution element). More elaborate treatments would 
incorporate a model for the small-scale gas clumping. For ex- 
ample, Barkana & Loeb (2004) model the density PDF based 



FIG. 2. — The transmitivity of a dynamic, infalling IGM surrounding a dark 
matter halo with total mass M tot = 5.2 X 1O U M0 at z = 3 (Left Panels) and 
Z = 5.2 (Right Panels). The horizontal axis shows the frequency offset from 
the line center, in units of the Doppler width. The upper/lower panels show 
the transmission in the absence/presence of a central ionizing source in the halo 
with Lj on = 5 X 10 44 ergs s , respectively. The solid-curve is the transmission 
from the center of our object, and is similar to the results found by Barkana 
& Loeb (2004). The dotted and dashed curves show the transmission at an 
angular impact parameter that corresponds to O.50 v j, and 1.00 v j,, respectively 
(9 v i r is the angular size of the virial radius of the dark matter halo). The figure 
shows the IGM transmission reaches a minimum redward of the true line center 
for all models. Furthermore, the IGM transmission is a function of impact 
parameter especially on the red side of the Lya line (—10 < x < 0), which may 
be important for a diffuse background source. 



on simulations. We find that our crude treatment yields compa- 
rable results for the transmission curves they obtained in their 
Figure 3 (especially when considering that their transmission 
will be reduced further near the line center in the absence of a 
quasar). 

The IGM transmission we obtain is presented in Figure 2 in 
four different cases. The transmittivity is shown as a function of 
frequency x for three different impact parameters, 9 = (solid- 
black curves), 9 = 0.59 vlr (red-dotted curves) and 9 = 1.0# v j r 
(blue-dashed curves). The upper panels show the transmission 
in the absence of any ionizing sources (apart from the uniform 
background), while the models shown in the lower panels in- 
clude a source with a total ionizing luminosity of 5 x 10 44 ergs 
s" 1 . We assumed a quasar-like power-law spectrum for the 
source of the form oc v~ x l . Barkana & Loeb (2004) performed 
similar calculations for the IGM transmission. The three main 
differences in the work presented in this paper are (i) our sim- 
pler treatment of the small-scale gas clumping, (ii) the possibil- 
ity that no bright ionizing quasar is present, and (iii) our spher- 
ically symmetric objects are spatially extended, and the IGM 
transmission depends on the impact parameter. The following 
can be inferred from Figure 2: 

• The IGM is most opaque at frequencies around the line 
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center. Depending on the exact model, the IGM's opac- 
ity is at its maximum slightly redward of the true line 
center. The reduced IGM transmission extends to fre- 
quencies corresponding to ~ v c ; lc , as already concluded 
by Barkana & Loeb (2003, 2004); 

• In the absence of an ionizing quasar, the transmission 
is reduced to ~ for a significant range of frequencies 
(whereas the transmission never drops to in Barkana 
& Loeb 2004, due to the presence of a bright quasar in 
their models). This range depends on the redshift, with 
attenuation decreasing at lower z, as shown by the differ- 
ence between the left and right panels. In the presence 
of the ionizing source, the attenuation is significantly 
reduced, as shown by comparing the upper and lower 
panels. 

• The IGM transmission is a function of impact param- 
eter. The frequency range with strong attenuation di- 
minishes at larger impact parameters, due to the smaller 
projected peculiar velocity of the IGM. 

2.4. Impact of a Dynamic IGM on the Results from Paper I. 

First, we show how the IGM modifies the spectrum emerg- 
ing from the fiducial model 1 from paper I. The left panel of 
Figure 3 shows the flux (in arbitrary units) as a function of fre- 
quency, x, for the intrinsic (blue-dotted line) and observable 
(black-solid line) cases. Since model 1. concerns the collapse 
of a neutral hydrogen cloud (i.e. no embedded ionizing source 
is present), the IGM transmittivity from Figure 2 for the 'no 
qso' case (upper panels) is taken. For this plot, the redshift of 
the Lya source was assumed to be z = 3. In the right panel, the 
flux is shown as a function of observable wavelength. The main 
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FIG. 3. — The impact of the IGM on the spectrum from our fiducial model 
1. The left panel shows the intrinsic (blue-dotted line) and observable (black- 
solid line) spectra. The IGM transmittivity from the upper left panel of Fig. 2 
is taken. In the right panel, the spectrum is shown as a function of observable 
wavelength. The figure shows that the IGM suppresses the blue peak and 
enhances the prominence of the dip between the red and blue peak. 

effects of the IGM here are to suppress the blue peak and to en- 
hance the prominence of the dip between the blue and red peak. 
We note that the dip remains at the same location. For the spec- 
tra emerging from the central models, the photons emerge with 
a larger average blue shift (with practically no photons redward 
of the line center), and for these models the main effect of the 
IGM is to suppress the total Lya flux. 

Second, we discuss how the IGM modifies the surface bright- 
ness profiles. In the absence of a quasar (or any other ionizing 



source), the IGM becomes more transparent with 6 for photons 
in the red peak. The result of this is that the surface brightness 
profiles of the reddest photons (x <, - 10) appears flatter than 
those calculated in paper I. This would strengthen the result 
that for neutral collapsing protogalaxies (i.e. in the absence of 
an ionizing source), the reddest Lya photons appear less cen- 
trally concentrated than the bluest Lya photons (see Figure 10 
in paper I). 

In the presence of a quasar (or any other ionizing source), 
the IGM also flattens the surface brightness profile of the red- 
dest photons (x SS - 10). However, the transmission of the 
IGM can decrease with increasing 9 for a range of frequen- 
cies, ^ x <, - 10, e.g. at x w -10 in the lower panels of Fig. 2. 
This is because at larger impact parameters, there can be a large 
effective column of absorbing material at projected velocity 
smaller than the maximum IGM infall speed. For photons in 
this frequency range, the IGM steepens the surface brightness 
profiles. 

The exact dependence of the IGM transmission on frequency 
and impact parameter in the presence (and absence) of an ion- 
izing source depends strongly on the exact assumed density and 
velocity field of the gas in the IGM. Furthermore, the IGM is 
unlikely to collapse in a spherically symmetric fashion. This 
introduces uncertainties to the angular dependence of the IGM 
transmission whose modeling is beyond the scope of this paper. 
Despite these uncertainies, the main impact of the IGM on the 
results obtained in paper I is clear: 

• The result that a steepening of the surface brightness 
profile towards bluer wavelengths within the Lya line 
is a useful diagnostic for gas infall, is not affected by 
the IGM. 

• The result that infall and outflow models can produce 
nearly identical spectra, may be affected by the IGM. 
The reason is that the IGM transmission is asymmetric 
around x = 0. Below (§ 3.1.2), we show that scattering 
in the IGM may cause emission lines that have identical 
profiles, but differ only by a small (~ several Doppler 
widths) frequency shift (caused by either inflow or out- 
flow), may have somewhat different observed values of 

(defined in paper I as the ratio of the number of 
photons in the blue and red peak). 

• All other results relating to (i) the net blueshift of the 
blue peak, (ii) FWHM of the blue peak, (iii) the location 
of the red peak (if any), (iv) the ratio [fi]/[7?] and its 
dependence on mass, gas infall velocity and redshift are 
not affected significantly by the IGM. 

With this knowledge at hand, and keeping in mind that a dy- 
namic IGM may produce an absorption feature slightly redward 
of the true Lya line center, we interprete several observations 
of known Lya emitters below. 

3. COMPARISON WITH OBSERVED LYMAN a SOURCES 

It will be useful to explicitly recall a result from paper I, in 
which we found the total observable Lya cooling radiation flux 
from an object of mass M tot , collapsing at redshift Zvu, to be 
approximately 
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where f a is the fraction of binding energy going into the Lya 
line. Other model assumptions included: 1) The gas density 
profile was assumed to trace that of the dark matter at large 
radii, but with a finite core. The emissivity profile was calcu- 
lated self-consistently by computing the rate of change in the 
gravitational binding energy at each radius; 2) The infall veloc- 
ity of the gas was taken to be of the form VbuikM = -v amp (r/ r v ; r ) a . 
The a dependence of /Lya is a fit that is valid for the range 
a <E [-0.5, 1] (see paper I for more details on our models). 

We caution that the strong break in the slope of the redshift 
dependence of the flux (eq. 4) at z = 4 should not be taken lit- 
erally. In reality, the variation of the this slope with redshift 
should be smoother. The main point of eq. (4) however, is that 
because the IGM becomes increasingly opaque towards higher 
redshifts at z > 4, cooling radiation of a fixed luminosity is more 
difficult to detect at z > 4 than at z ^ 4. Quantitatively, the total 
detectable Lya cooling flux does not scale with (which can 
be approximated as oc (1 + z)~ 2 - 75 ), but as oc< e~ T > dj 2 (which 
can be approximated as oc (1 +z)~ n ) for z > 4.0. 1 

In paper I, we studied a range of models with different masses, 
redshifts, and emissivity, density, and velocity profiles. Table 2 
summarizes the parameters of two of paper I's fiducial models 
and two models used in § 3.2. For each model the table contains 
the total mass, redshift, and circular velocity of the host halo, as 
well as the parameters describing the infall velocity field. For 
each model, all gas continously emits all its gravitational bind- 
ing energy as Lya as it navigates down the dark matter potential 
well (in the terminology of paper I, these are 'extended' models 
with f a = 1.0). 



Table 2 
Model Parameters 



Model # 


M tot 


Zv'u 


Vcirc 


Vamp 


a 


fc 


(Ext.) 


1O U M 




km s~' 


Vcirc 






1 


5.2 


3.0 


181 


1.0 


1.0 


1.0 


2 


5.2 


3.0 


181 


1.0 


-0.5 


1.0 


Ml 


15 


3.0 


257 


1.0 


1.0 


1.0 


M2 


15 


3.0 


257 


1.0 


-0.5 


1.0 



3.1. Comparison with Steidel etal's LABs. 

We start our comparison with the data presented in Steidel 
et al. (2000), who found two large, diffuse, and relatively bright 
Lya emitters ('Lya blobs' or 'LABs') in a deep narrow band 
search of a field known to contain a significant overdensity of 
Lyman Break galaxies. The total observed flux is ~ 10~ 15 ergs 
s" 1 cm" 2 for each blob, and both have a physical size exceed- 
ing 140 kpc. The excitation or ionization mechanism for either 
blob is poorly understood. One blob, referred to as 'LABI', 
has a bright sub-mm source as a counterpart (with some uncer- 
tainty in the identification, given the 1 5" half power beam width; 
Chapman et al. 2004). The physical interpretation is that a dust 



'This conclusion does not account for the proximity effect near bright QSOs, 
which may allow significantly more flux to be transmitted (see Paper I for more 
details). 



enshrouded AGN, which is undetected at other frequencies, is 
buried within the Lya halo. It has been speculated that a jet-like 
structure that is well outside of our line of sight induces star for- 
mation, which, combined with the ionizing radiation from the 
AGN, photo-ionizes the surrounding hydrogen cloud (Chapman 
et al. 2004). The other blob, 'LAB2', has no sub-mm bright 
counterpart, but Basu-Zych & Scharf (2004) have found an as- 
sociated hard X-ray source in Chandra archival data, with a to- 
tal luminosity comparable to that of the Lya (L <~ 10 44 ergs s" 1 ). 
Basu-Zych & Scharf (2004) have speculated that both LABI 
and LAB2 may have an embedded AGN, which goes through 
periodic episodes of activity. An old population of relativis- 
tic electrons, originating near the supermassive black hole that 
powers the AGN, can upscatter CMB photons into the X-ray 
band, which may then illuminate the hydrogen gas and explain 
at least part of the Lya flux in these blobs. This scenario is 
based on interpretations of recent findings of extended X-ray 
sources at redshifts as high as z = 3.8 (e.g. Scharf et al. 2003). 

The possibility that cooling radiation contributes significantly 
to Lya flux of the LABs remains. First, we notice that to pro- 
duce a flux of 10~ 15 ergs s" 1 cm -2 from cooling radiation alone, 

requires v circ ~ 300 km s _1 orM tot ~ 2-3 x 10 12 /a 3/5 M Q . The 
virial radius of an object of this mass is ~ 100 fa ^ kpc, corre- 
sponding to ~ 13 fa I 5 ", in good agreement with the observed 
size of the blobs. Given that the image can be observed over 
a large region on the sky suggests that the surface brightness 
profile is not declining rapidly with radius, which argue against 
our "central" model. More quantitatively, Figure 9. of paper I 
shows that d\nF I d\nQ <^ -2 and <; - 1 for the central and ex- 
tended models, respectively for M = 4 x 1O 12 M . The plots in 
Steidel et al. (2000) lack contour levels, but especially the sur- 
face brightness of fainter emission does not appear to rapidly 
decrease with radius. 

For extended cooling radiation, the Lya spectrum should be 
double peaked (provided the red peak is not swallowed by the 
blue, but for v amp ~ v c ; rc ~ 300 km s" 1 , this is not expected ac- 
cording to our findings in paper I). Steidel et al. (2000) present 
spectra taken along slits across both blobs. More recently, how- 
ever, Bower et al. (2004) and Wilman et al. (2005) have used the 
SAURON Integral Field Spectrograph to observe Steidel et al. 
(2000)'s LABI and LAB2, respectively. Since these spectra 
cover the entire two-dimensional region of Lya emission, we 
use these more recent spectra for comparisons to our models. 

3.1.1. Comparison with Spectra of LABI 

The spectra presented by Bower et al. (2004) have a complex 
structure. If LABI is emitting cooling radiation, surrounded 
by an extremely optically thick (t <; 10 6 ) collapsing hydrogen 
cloud, the Lya radiation is expected to emerge with a net blue 
shift. It is not possible to tell from Bower et al. (2004)'s data 
whether this is the case, since, in the absence of other lines, 
the systemic redshift of the LAB is not known. There is no evi- 
dence for a fainter redder peak. In Steidel et al. (2000)'s original 
data, there is evidence that the main bluer peak is accompanied 
by a fainter redder peak, with ^ 1. Unfortunately, it is 

not possible to calculate the exact value of from Steidel 

et al. (2000)'s paper, since the contour levels of their spectra 
are not shown. In our fiducial model, but with the mass of the 
halo increased to M tot = 3 x 10 12 M Q , as required by the angular 
size and total flux from LABI, we expect ^100 after 

including the effect of the IGM opacity (see Figure 8 in paper 
I for how the [B]/[R] ratio scales with various parameters). A 
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secure detection of a double-peaked, asymmetric profile, with 
a stronger blue peak, [.B]/[/?] 3> 1, would be strong evidence 
for an optically thick, collapsing gas cloud. Bower et al. (2004) 
note that the line width appears to increase with radius, a feature 
that does not arise in models of neutral collapsing gas clouds 
presented in paper I. However, it does arise in models in which 
an ionized source is embedded within the gas cloud (see Fig. 
13 of paper I). Indeed, the presence of a sub-mm source sug- 
gests that an ionizing source resides within LAB 1 . The ionizing 
source can not, however, fully ionize the gas (Fig. 13 in paper 
I), without violating [.B]/[/J] ^ 1 as observed by Steidel et al. 
(2000). 

The origin of the Lya emission in our model is then attributed 
to fluorescent radiation emitted in an HII region created by the 
central ionizing source. This HII region is surrounded by a 
neutral collapsing hydrogen envelope. Cooling radiation, gen- 
erated over a spatially extended region within the neutral hy- 
drogen envelope, may provide an additional contribution to the 
Lya luminosity. Evidence for this scenario would be (i) a flat 
surface brightness profile, (ii) the strong asymmetry of the double- 
peaked profile ([B]/[R] > 1) and (iii) the increase of the FWHM 
of the line with radius. This interpretation is consistent with that 
given by Bower et al. (2004). 

We note that Ohyama et al. (2003) have also obtained deep 
spectroscopy in a slit across LABI, and showed that the spec- 
trum of the brightest knot in particular, is more complex than 
can be explained with our models. The spectrum of the cen- 
tral knot can be be fitted with three Gaussians, each separated 
by ~ 2000 km s -1 . The other regions are consistent with our 
model above. 



3.1.2. Comparison with Spectra ofLAB2 

LAB2 has a complicated morphology, with six distinct sub- 
blobs of resolved Lya emission. The spectra of the sub-blobs 
all exhibit a double-peaked frequency structure. Wilman et al. 
(2005) show that the spectra can be explained with a very sim- 
ple model. The emission line is assumed to be a Gaussian with 
a FWHM of ~ 1000 km s" 1 , with its peak frequency varying 
by ~ 290 km s" 1 over the entire LAB. A foreground slab of 
neutral gas with Nhi ~ 10 19 cm" 2 leaves an absorption feature 
which falls at the same wavelength in all six spectra. Wilman 
et al. (2005) attribute this absorption feature to a foreground 
shell of neutral hydrogen gas, swept up by a galaxy-sized out- 
flow. Indeed, in four out of six sub-blobs (except for B and F), 
the absorption line lies blueward of the peak of their Gaussian 
emission line. We show below, however, that these data can 
also be explained with an equally simple in-fall model. 

In our model, we take the same intrinsic emission profile 
of the Lya emission as in Wilman et al. (2005), and likewise, 
we allow ourselves the freedom to shift the peak frequency of 
the intrinsic emission to the blue by up to ~ 300 km s" 1 (see 
the shift assumed in each panel). We process this emission 
line through the model of the dynamic IGM as discussed in 
§ 2. We calculated the IGM transmission curves expected for a 
source with M tot = 10 12 M Q and added an ionizing source with 
Li on = 2 x 10 44 ergs s" 1 to increase the transmission near the 
line center (Note that we assumed 8 = in eq. 2 for all panels, 
as discussed in § 2). The results are shown in Figure 4, which 
can be compared directly to Figure 3 in Wilman et al. (2005). 
Our model provides an equally good fit to the data. Further 
improvement of the fits is possible, if we allow ourselves the 



freedom to vary the IGM model (and the exact shape of the 
emission line) from location to location on the sky, rather than 
imposing a single IGM infall model as we have done. 

To understand why our model provides a good fit is through 
the IGM transmission curve shown in lower left panel of Fig- 
ure 2 2 . Panel D in Figure 4 is obtained by multiplying a Gaus- 
sian emission line that is centered on the line center (x = 0), by 
this IGM transmission. Shifting the Gaussian back and forth 
in frequency and multiplying these emission lines by this trans- 
mission curve yields the results shown in Figure 4. Each panel 
contains the required velocity shift in its top right corner. A 
negative velocity corresponds to a blue shift of the line. 
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FIG. 4. — A reproduction of the observed Lya spectra of Steidel et al's 
LAB2 by Wilman et al. (2005). The spectra in the different panels are dis- 
played in the same way as originally in Wilman et al. (2005), who succesfully 
reproduced their observations with a foreground shell of neutral hydrogen gas, 
swept up by a galaxy-sized outflow. This figure demonstrates that our infall 
model reproduces the observed spectra of LAB2 equally well. Our model as- 
sumes simple infall for the IGM, as outlined in § 2. The red-solid lines are 
the spectra after scattering by the IGM. The solid-black lines are the intrinsic 
spectra. The spectra in panels B and F, with enhanced blue peaks, are consis- 
tent with cooling radiation that is propagating outwards through the collapsing, 
neutral hydrogen gas. 



First we point out that all panels, except for panel D, require 
an intrinsic blue shift of the emission line. This requires the gas 
within the virial radius to be collapsing. Second, the spectra in 
panels B and F, which are the two farthest out-lying sub-blobs 
in the two-dimensional image of LAB2, have enhanced blue 
peaks, requiring a significant blueshift in the intrinsic emission, 
by ~ 3-400 km s _1 . Within our model, this blue-shift would 
naturally arise if the gas far out at the locations of sub-blobs B 

2 We point out that in Figure 4 the horizontal axis displays the observable Lya 
wavelength (blue and red lies on the left and right, respectively), while the 
horizontal axis of Figure 2 displays a frequency variable (and, consequently, 
blue and red lies on the right and left, respectively). 
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and F is collapsing and optically thick (r Q > 10 4 ) to Lya pho- 
tons. The smaller intrinsic required blueshift in panels A, C 
and E can be accounted for if the Lya photons were propagat- 
ing outwards through gas that is photoionized by an embedded 
ionizing source (see Figure 1 1 in paper I). Wilman et al. (2005) 
show in their Figure 2 that regions B and F are also signifi- 
cantly fainter than the other sub-blobs. This would again be 
expected the Lya radiation emerging from these regions were 
attributed to cooling radiation (as opposed to photoionization 
for the other sub-blobs, which reside closer to any embedded 
ionizing source). Our physical picture associated with the col- 
lapse is that the gas is collapsing as a whole onto their mutual 
center of mass. Panels B & F show Lya cooling radiation, gen- 
erated over a spatially extended region. The Lya emission from 
the blobs in the other panels consists of processed ionizing radi- 
ation from embedded ionizing sources. The IGM surrounding 
LAB2 is dragged along in the overall cloud collapse. 

Based on the data on the two Lya blobs presented in Steidel 
et al. (2000), we conclude that the spectrum and surface bright- 
ness distribution of LAB 1 are consistent with a model in which 
gas is collapsing in a DM halo with M tot = a few xlO 12 M 
while emitting fluorescent Lya radiation in response to pho- 
toionization by (an) embedded ionizing source(s) (except the 
brightest knot, which is likely powered by an outflow Ohyama 
et al. 2003). We found that the spectra of LAB2 can be ex- 
plained if the intrinsic Lya emission profile is a single Gaus- 
sian, blue-shifted by up to several hundred km/s due the col- 
lapse of optically thick gas within a massive halo, and subse- 
quent resonant scattering on larger scales within a perturbed 
IGM. 

3.2. Comparison with Matsuda et al.'s LABs 

More recently, Matsuda (2004) published a catalogue of an 
additional 33 LABs, discovered in a deep, narrow-band, wide 
field (31' x 23') search with Subaru, centered on the same 
field as studied by Steidel et al. (2000). One third of their 
blobs are not associated with bright UV sources, and thus pro- 
vide good candidates for objects dominated by cooling radia- 
tion. Mori et al. (2004) have used ultra-high resolution simula- 
tions to model a forming galaxy undergoing multiple epochs of 
supernova explosions and found their models could reproduce 
the typical Lya luminosities observed from these LABs. Their 
models are particularly successful in reproducing the 'bubbly' 
appearance of several LABs. Matsuda (2004) suggest that five 
specific Lya blobs may be cooling objects, based on the dif- 
fuse appearance of the images. A diffuse appearance, may not 
be a necessary requirement for cooling radiation. The physi- 
cal picture associated with extended cooling is that clumps of 
cold gas are flowing inwards, while emitting Lya in response 
to various sources of heating (i.e. gravitational shock heating 
and possibly photoionization heating by the thermal emission 
from the hot, ambient gas). Therefore, some dumpiness in the 
surface brightness distribution may be expected. If we drop the 
requirement that the Lya image is diffuse, then the number of 
candidates for cooling objects is increased to 12. 

The detected fluxes for the majority of the LABs lie in the 
range 7 x 10~ 17 -2 x 10~ 16 ergs s" 1 cm" 2 , which correspond to 
cooling radiation from objects in the mass range M tot ~ 1 — 2 

x lO 12 / a 3 ^ 5 M0, for extended models with (v amp ,a) = (v c ; rc ,l) 
(see eq. 4). The corresponding range in virial radii is ~ 80- 
lOO/a'^ 5 kpc, which translates to an angle on the sky of ~ 10— 
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FIG. 5. — The figure shows the surface brightness profile of a model that 
has the observed gross properties of the majority of the blobs found by Mat- 
suda (2004) (black solid line, see text). If this model is correct, then deeper 
observations will reveal these blobs to be spatially more extended. The red 
dotted line represents the same model with a changed from —0.5 to 1.0. This 
object would escape detection since its surface brightness is everywhere below 
Matsuda (2004)'s detection limit, indicated by the thick-dotted horizontal line. 

12/ Q 1//5 ". Converting this to an area yields values which exceed 
the observed Lya isophotal area by a factor of ~ 10. This is not 
surprising since the calculated surface brightness distribution 
may drop below the detection threshold in the outer parts. This 
is illustrated in Figure 5, in which we plot the surface brightness 
profile for model Ml as the red-dotted line. This is equivalent 
to our fiducial model, but with M tot increased to 1 .5 x 10 12 Mq . 
The total flux on earth for this model is predicted to be 1 .6 x 
10~ 16 ergs s" 1 cm" 2 . 

From Figure 5, however, it can be seen that this radiation 
is too diffuse to be detected. The total flux is spread over an 
area of ~ 300 arcsec 2 , which brings the average surface bright- 
ness below the detection threshold of 3 x 10~ 18 ergs s" 1 cm ~ 2 
arcsec" 2 . The flux is more centrally concentrated for models 
with lower a (see paper I). Plotted as the black-solid line, is 
the surface brightness profile of model M2, which is identi- 
cal to model Ml, except with a = -0.5. For this model, the 
flux from the inner 3" exceeds the detection limit. The total 
flux from this detectable region is 1.5 x 10~ 16 ergs s" 1 cm -2 , 
which is spread over an area of ~ 30 arcsec 2 . This angular ex- 
tent is consistent with the observed area of the majority of Lya 
blobs, 20-40 arcsec 2 . The average surface brightness of the 
detectable region of model M2 is ~ 5.2 x 10~ 18 ergs s" 1 cm ~ 2 
arcsec" 2 , which corresponds to mab ~ 27.1 mag arcsec" 2 . Mat- 
suda (2004) focussed on imaging and did not present spectra in 
their paper, so further comparison is at this stage not possible. If 
this model is correct, then deeper observations will reveal these 
sources to be spatially more extended. We showed in Figure 10 
of paper I that for extended models with a = -0.5, the bluest 
15% of the Lya photons would appear more centrally concen- 
trated than the reddest 15%. Comparing Figure 10 from paper I 
and Figure 5 suggests that the surface brightness of the reddest 
15% of the photons may have been below the detection limit, 
implying that deeper observations may reveal the red peak in 
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the spectrum (if it is not already in the data). 

3.3. Spectroscopic Detections 

Our comparisons have so far focused on extended LABs, 
which are a rarity among the known Lya emitters. As men- 
tioned in § 1, narrow band searches on Kitt Peak (LALA, e.g. 
Rhoads et al. 2000; Rhoads & Malhotra 2001; Rhoads et al. 
2004), on Keck (e.g. Hu et al. 1998) and on Subaru (e.g. Ko- 
daira 2003; Hu et al. 2004; Ouchi 2005), have been very suc- 
cessful in finding galaxies at redshifts z > 4.5. The narrow 
band surveys, combined with deep broad band images, pro- 
vide candidates based on selection criteria involving the col- 
ors, equivalent widths, and flux thresholds. These candidates 
need spectroscopic confirmation, using large telescopes such as 
the Very Large Telescope (VLT) and Keck. Roughly 30% of 
the candidates are confirmed to be genuine high redshift Lya 
emitters (e.g. Hu et al. 2004). Contrary to the lower redshift 
LABs, these Lya emitters thus all have observed spectra. The 
known Lya emitters typically have total fluxes in the range 
5 x 10~ I8 -5 x 10~ 17 ergs s" 1 cm" 2 , which is fainter by at least 
an order of magnitude than the LABs discussed above. 

If some of these Lya emitters are cooling objects, then deeper 
observations will reveal fainter, more diffuse emission associ- 
ated with them (unless the emission is strongly centrally con- 
centrated, in which case the surface brightness profiles drops 
rapidly, see Fig. 1 1 of paper I). To produce a detectable flux 
of 10~ 17 ergs s" 1 cm" 2 from a z = 5.7 source requires M tot ~ 
1O 12 M (eq. 4). A dark matter halo of this mass that col- 
lapses at z = 5.7 is associated with a ~ 3a fluctuation in the pri- 
mordial density field. According to the Press-Schechter mass 
function, the number density of halos with M > 10 12 M©, is 
1.3 x 10" 5 Mpc" 3 at z = 5.7, which is factor of ~ 40 lower 
than the implied number density of Lya emitters at this red- 
shift (e.g. Ouchi 2005). Additionally, the typical observed flux 
from a Lya emitter of 10~ 17 ergs s" 1 cm" 2 comes from a small 
region, ~ 1 arcsec 2 . This implies that either all cooling radi- 
ation from this M tot = 10 12 M Q object emerges from this small 
area, which requires the surface brightness profile to be very 
steep, and thus favor the central models (§1). Alternatively, we 
could be detecting the brightest central region a spatially ex- 
tended Lya halo (as in Fig. 5). This, however, would require an 
even more massive, and therefore even rarer, halo to produce 
the measured Lya flux levels. The above suggests that the ma- 
jority of the known high redshift Lya emitters are not powered 
by cooling radiation alone. Indeed, the continuum detection of 
many high redshift Lya emitters, which show that the equiva- 
lent width of the Lya line is consistent with ionization by young 
stars forming with a usual IMF. 

An object of particular interest to us is described by Bunker 
et al. (2003b), who discovered a z = 5.78 galaxy with a total 
Lya flux of 2 x 10~ 17 ergs s" 1 cm" 2 in the Chandra Deep Field- 
South. Using the DEIMOS spectrograph on Keck, they obtained 
a spectrum, which looks remarkably like the spectrum for a 'ex- 
tended' model: a small red peak is separated from a more pro- 
nounced bluer peak by ~ 400 km s" 1 . If powered by cooling 
radiation, a total flux of 2 x 10~ 17 ergs s" 1 cm" 2 would require 
gas cooling in a halo with a mass of M ~ 10 12 M Q . If the gas in 
this object were continuously cooling, this would produce the 
maximum of the red peak to lie > 300 km s" 1 from the line 
center. The problem with this interpretation however, is that 
the Lya emission is extremely compact, with an angular size 
less than ~ 0.1". Since the emission is detected at the 20 -a 
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FIG. 6. — This figure shows a the spectrum of a model that is consistent 
with recent observations by of a z = 5.78 Lya emitter in the Chandra Deep 
Field— South by Bunker et al. (2003b). An infalling IGM, as outlined in § 2, 
can modify the intrinsic Gaussian emission line (indicated by the black-dotted 
line) to have a dip in the observable spectrum (indicated by the red-solid line) 
on the red side of the line. 



level, this implies that the surface brightness falls of by a factor 
of <; 20 over ~ 0.1". Because extended cooling models have 
rather flat surface brightness distributions (with dlnF I d\n0 
~ -1 to -0.5, see paper I), they are not consistent with the ob- 
served compactness of the emission. A similar conclusion ap- 
plies to the faint, strongly gravitationally lensed Lya emitter 
discovered recently behind a galaxy cluster (Ellis et al. 2001), 
which remains a point-source despite the large increase in the 
effective spatial resolution afforded by the presence of strong 
lensing. 

An alternative interpretation of this spectrum combined with 
the compactness of the emission is that the Lya emission is due 
to stars in an already formed galaxy, as Bunker et al. (2003b) 
argue, that fully photoionizes the gas in the galaxy. The emis- 
sion line is therefore centered on the true line center and its 
FWHM set by the gas' velocity dispersion. If we process this 
Lya through an infalling IGM, as described in § 2, we find we 
can reproduce the gross features of the observed spectrum, as 
shown in Figure 6 by the solid-red line. The intrinsic line pro- 
file is given by the black-dotted line. In this model, the IGM 
scatters ~ 95% of the Lya photons out of the line of sight. It 
is generally assumed that the IGM only allows photons from 
high-z Lya emitters to be detected on the red half of the Lya 
emission line, implying ~ 50% is scattered out of our line of 
sight. If our model is correct then the intrinsic luminosity may 
of this Lya emitter may have been underestimated by a fac- 
tor of ~ 10. The intrinsic Lya equivalent width would then be 
increased from 20 to ~ 200 A. This could have important con- 
sequences: if the intrinsic Lya luminosity has been systemati- 
cally underestimated for high-z Lya emitters, then the escape 
fraction of Lya photons would have to be unusually high (for a 
constant star formation rate). This may put interesting limits on 
the dust abundance and/or distribution in these galaxies. Alter- 
natively, it could imply that estimates of the total intrinsic Lya 
luminosity must be revised, requiring a ~ 10 times larger star- 
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formation rate. The latter would enhance the implied contribu- 
tion of high redshift Lya emitters to the ionizing background. 
Either way, these issues deserve more attention and will be ad- 
dressed in more detail elsewhere (Dijkstra et al. 2006). We em- 
phasize that the only requirements for these issues to be relevant 
for a Lya emitter are that i) its surrounding IGM is collapsing, 
and ii) its intrinsic emission line is not redshifted significantly 
(Vshift Vcirc, but note this can easily be achieved in outflows). 

4. DISCUSSION 

As argued above, observational evidence exists for gas in- 
fall around both pointlike (§ 3.3) and extended (§ 3.1.2) Lya 
emitters. Below (§ 4.1) we compare our model for Steidel 
et al. (2000)'s LAB2 with a recently proposed outflow model 
by Wilman et al. (2005). We discuss why two such different 
mechanisms can reproduce the observed spectra well. We argue 
in § 4.2 that the notion that infall models may produce similar 
Lya spectral features as outflow models may be more generic 3 
and discuss its cosmological implications. In § 4.3 we discuss 
how additional imaging in Ha may help to resolve this issue. 

4. 1 . Infall vs. Outflow in Steidel et al's LAB2 

It may come as a surprise that two completely different mech- 
anisms (infall vs. outflow) may reproduce the observed spectra 
of Steidel et al. (2000)'s LAB2. To produce a spectrum with 
< 1 using a single absorbing shell of material, the shell 
must absorb on the blue side of the emission line, and must 
therefore be outflowing. As mentioned above, an infall model 
naturally produces < 1 (shown in Fig. 2). Note that the 

for the exact same intrinsic emission line and the same 
ratio, the dip in the spectra for the outflow and infall model lies 
blueward and redward of the true line center, respectively. This 
difference is in practice difficult to see, especially since the in- 
trinsic emission line is unlikely to be Gaussian (see Fig. 3 and 
paper I, in which we calculated the exact shape of the emission 
line emerging from models of both neutral and photoionized 
collapsing gas clouds). We caution that the isotropically in- 
falling IGM will not change the total flux toward Earth. How- 
ever, the scattering "IGM shell" is geometrically quite thick (up 
to ~ 5 virial radii ), so any photons scattered back into our line 
of sight will be spread over a large solid angle (and also spread 
in frequency). This "super-fuzz" is too diffuse to modify the 
simple absorption dip due to scattering out of the line of sight. 

An advantage of our explanation is that the intrinsic blue shift 
of the emission line, in all panels, is consistent with overall gas 
collapse. The amount of blueshift of each line is determined by 
the brightness of the photoionized source. Indeed, the brightest 
Lya blobs have the smallest net blue shift. 

A second advantage of the above explanation of the data 
is that it may require less fine tuning than the explanation by 
the swept-up shell, given by Wilman et al. (2005). The peak 
frequency of the absorption dip varies by less than 60 km s -1 
among the six different sub-blobs covering LAB2. This small 
variation would arise naturally if the entire shell has virtually 
stalled and come to rest with respect to the galaxy's systemic 
velocity. Because a superwind will spend most of its time close 
to the zero velocity stalling radius it is most likely observed at 
this position. However, for the shell to stall and preserve a sim- 
ilar column density over large length scales (a few 100 kpc) and 
long times (10 8 yrs), requires the ambient pressure of the IGM 

3 Note the similarity of the Lya spectrum produced by outflow and infall mod- 
els is mainly due to the IGM and is therefore unrelated to our result (iv) in § 1 



not to exhibit pressure gradients on these scales. In our model, 
the IGM infall velocity is <~ 300 km s" 1 and may vary by ~ 20% 
over the entire blob and may be a less stringent requirement. 

In either model, the difference between LAB 1 and LAB2 re- 
mains a puzzle. If the shape of the observed spectrum in LAB2 
is indeed due to IGM infall, then this IGM infall is also ex- 
pected around LAB 1 . However, no evidence for this is seen in 
the spectra of in LABI. In our infall model, this would require 
that the Lya flux emerges across LABI with a larger intrin- 
sic blueshift of ~ 500 km s" 1 . As we showed in paper I, this 
can be easily reached when the gas in LAB 1 is optically thick 
(T a ^ 10 4 ) to the Lya photons. 

4.2. Infall vs. Outflow in all Lya Emitters 

Ample spectroscopic evidence exist for outflows of gas around 
Lyman Break Galaxies (LBGs, e.g. Adelbergeret al. 2003, and 
references therein). It is not yet securely established to what 
level these superwinds surround the galaxy (i.e. what solid an- 
gle around the galaxy is swept up by the wind). The physical 
extent of the superwind is set by the star formation history of 
its host galaxy, and determines to what degree it can prevent 
gas from the IGM to continue to cool and collapse onto the 
galaxy. Understanding superwinds is clearly a requirement for 
understanding galaxy formation and evolution. Wilman et al. 
(2005) have argued that their recent observations demonstrate 
that superwinds occur on galactic scales around Steidel et al. 
(2000)'s LAB2. However, we showed in § 3.1 above that these 
observations are also consistent with an inflowing IGM. 

One may argue that the number of LBGs with known out- 
flows greatly outnumber the number of known inflows. How- 
ever, this does not necessarily imply that galaxies with outflows 
are more common than those with significant gas infall. At 
high redshifts, z > 5, candidates for galaxies containing a su- 
perwind are found based on their Lya properties (see Taniguchi 
et al. 2003, for a review). As in the case of LAB2, some of the 
Lya spectra of candidate superwind galaxies (e.g. Dawson et al. 
2002) could in fact also be explained by infall of the IGM (Di- 
jkstra et al, 2006b, in prep) onto the Lya emitter. Additionally, 
for a given intrinsic Lya luminosity, there is observational bias 
to detect an outflow. The intrinsically redshifted Lya emission 
from outflows is easier to detect than the blueshifted Lya from 
infall, because the latter will be subject to much more severe at- 
tenuation in the IGM. Furthermore, outflows can be more com- 
pact than large scale infall, increasing the surface brightness. 
These observational biases must be accounted for when calcu- 
lating the fraction of Lya emitters associated with outflows and 
infall. With ongoing and future Lya surveys, the sample of high 
redshift Lya emitters steadily grows, allowing a better determi- 
nation of this number, and addressing various basic questions 
related to galaxy formation: (i) how common is Lya cooling 
radiation among collapsing protogalaxies? (ii) how common 
is fluorescent Lya emission due to embedded ionizing sources 
(such as a central quasar or young stars that formed during the 
collapse) among protogalaxies? iii) how common are super- 
winds in young galaxies at high redshifts? 

4.3. Ha Emission from Protogalaxies 

Discussion of emission line radiation from protogalaxies have 
focused typically on Lya radiation, but as argued by Oh (1999), 
the Ha line may be a cleaner probe of the ionizing emissivity. 
The production rate of Ha photons is 0.45 per recombination, 
which is comparable to that of Lya. Because hydrogen atoms 



10 



spend such a short time in their n = 2 state, the Ha photons do 
not undergo resonant scattering, and escape from protogalax- 
ies in a single flight. Because an Ha photon is 5/27 times as 
energetic as Lya, the total intrinsic Ha flux is lower. How- 
ever, since Ha photons do not suffer any attenuation in the IGM 
(Oh 1999), the detectable Ha flux can exceed that of Lya, es- 
pecially for neutral collapsing gas clouds at higher redshifts. 
A simultaneous detection of both Ha and Lya would allow a 
separation of the intrinsic emissivity from radiative transfer ef- 
fects. For example, a relative red or blueshift of the Lya line 
compared to the Ha line is strongly indicative for an outflow or 
infall, respectively and separates infall from outflows. 

Assuming that the intrinsic Ha flux is lower than the Lya 
flux by factor of Tua ~ 27/5 = 5.4, we estimate the total de- 
tectable Ha flux emerging from a neutral collapsing gas cloud 
(eq. 4) to be: 

ha / 2-a \'- 2 / 5 w v circ \ 5 / 5 \ 2 - 75 ergs 

10" 18 ~ ■ U.75/ \F H J Uhlans" 1 / Vl+zJ s cm 2 ' 

(5) 

where we used v amp = v c i rc and that f Ha = L Ha /[Amij}. A Fluxes 
at this level will be detectable by the James Webb Space Tele- 
scope (JWST). The Near Infrared Spectrograph (NIRspec) on 
JWST can detect a ~ 100 nJy point source at the lOer level in a 
10 5 s exposure time 5 . The flux calculated in eq. (5) is spread 
over a frequency range of Av ~ 2v c i rc fH Q (l +z)/c. This trans- 
lates to a flux density of 

hay _ / 2-a \ L2 / 5 W v circ \ 2 / 5 \ L75 

1400 \ 1.75/ 1.7^/1116 km s" 1 / Vl+J " Y ' 

(6) 

which is spread over an area 7r# 2 ir , where f? V ir = r^/dAiz), in 
which d A is the angular diameter distance. Given that a point 
source in NIRspec subtends only ^0.01 arcsec 2 on the sky, 
the total noise over the area over which the Ha is emitted is 
[(7r6» 2 r )/(0.01 arcsec 2 )] 05 times higher. This translates to the 
following signal to noise ratio in a 10 5 s exposure: 

— ~ in x ( ^ Ha ' v \ ( °' 01 arcsec2 ^| 1/2 ~ 
N~ UOOnJy/V n9 2 vh / 

which is a very weak signal. However, the Ha emission is 
not uniformly emerging from the area enclosed by v ; r , but is 
(much) higher in the central region of the image, which can 
boost the expected flux densities. For model 1. and 2. (Ta- 
ble 2) for example, 50% of the flux comes from 9 <^ 0.27# v j r 
and 9 <> 0.18# v ir, respectively. This would boost the S/N for 
model 1. and 2. by a factor of <~ 2 and 2.5, respectively. It 
should be noted that with the current proposed spectral resolu- 
tion of R = 1000, the Ha line would not or barely be resolved. 
Because of increased detector noise at A > 5 /xm, the S/N ratio 
drops by almost an order of magnitude for Ha at z > 7. Oh 
et al. (2001) argued that because the helium equivalent of Ha, 
which has a transition at A = 1640 A, does not suffer from the 
increased detector noise for z ^ 32, it may be more easily de- 
tectable. This would only apply to our self-ionized cases, and 
under the restriction that the spectrum of the ionizing sources is 
sufficiently hard (3 ~ 1 (see Oh et al. 2001). 

4 We approximated the luminosity distance by Jl(z)=3.6 x 10 4 [(1 +z)/5] 11//8 
Mpc as in § 3. 

5 See http://www.stsci.edu/jwst/science/sensitivity. 



5. SUMMARY AND CONCLUSIONS 

In paper I, we calculated the properties of the Lya radiation 
emerging from collapsing protogalaxies. In this paper we pro- 
vided simple calculations of the mean transmission of the IGM 
around such objects to facilitate a comparison with observa- 
tions. We assumed that the IGM surrounding our objects to be 
flowing in, following the recent prescription by Barkana (§ 2). 
We found the main results from paper I not to be affected sig- 
nificantly. 

In § 3 we compare our models with the observations of known 
Lya blobs presented by Steidel et al. (2000) in § 3.1, and Mat- 
suda (2004) in § 3.2. The observed properties of those LABs 
without significant continuum in the larger sample of blobs pre- 
sented by Matsuda (2004), such as the surface brightness and 
angular size, are consistent with model in which gas is cool- 
ing and condensing in dark matter halos with M to t ~ 1 — 2 x 
10 12 M Q . According to this model, deeper observations should 
reveal that these source are spatially extended by a factor of 
several beyond their currently imaged sizes (see Fig. 5). With 
the current data it may already be possible to distinguish that 
the bluest 15% of the Lya photons would appear more cen- 
trally concentrated than the reddest 15%, which we propose as 
a diagnostic of optically thick gas infall. 

Apart from the brightest central knot, which Ohyama et al. 
(2003) have shown to be likely powered by an outflow, the ob- 
served spectrum and surface brightness distribution of Steidel 
et al. (2000)'s LABI are consistent with a model in which gas 
is collapsing in a DM halo with M tot = a few x 10 12 M Q . An em- 
bedded ionizing source photoionizes the central part of the col- 
lapsing cloud, in which fluorescent Lya is emitted. Cooling ra- 
diation, emitted in the surrounding neutral hydrogen envelope, 
may provide an additional contribution to the Lya luminosity. 
We found that the spectra of LAB2, as presented recently by 
Wilman et al. (2005), can be explained if the intrinsic Lya emis- 
sion profile is a single Gaussian, blue-shifted by up to several 
hundred km/s due the collapse of optically thick gas within a 
massive halo, combined with subsequent resonant scattering on 
larger scales within a perturbed, infalling IGM. 

We found similar evidence for infall of the IGM onto the 
z = 5.78 Lya emitter observed by Bunker et al. (2003b) in its 
spectrum (§ 3.3). We argue that IGM infall around high redshift 
Lya emitters could be a more common phenomenon than previ- 
ously believed, as infall models may produce spectra similar to 
those of outflow models. If this is correct, then the intrinsic Lya 
luminosities, and derived quantities, such as the star-formation 
rate of high redshift Lya emitters, may have been underesti- 
mated significantly (§ 4.2). 

Because of the increasing opacity of the IGM with redshift, 
we expect the detectable Lya fluxes from high redshift (z > 4.0) 
collapsing protogalaxies to be very weak, unless the intrinsic 
Lya luminosity is very high. We concluded that the majority of 
the observed high redshift Lya emitters (§ 3.3) are most likely 
not associated with collapsing protogalaxies. At high redshift 
(z > 4.0) the detectable Ha flux emerging from neutral col- 
lapsing gas clouds may dominate that of Lya (§ 4.3). In cases 
in which both the Lya and Ha are detected simultaneously, it 
would be possible to separate the intrinsic emissivity from ra- 
diative transfer effects. 

The work in this paper emphasizes the importance of the Lya 
line as a probe of the earliest stages of galaxy formation. In par- 
ticular, with the growing sample of high redshift Lya emitters 
(both point like and extended), and with deeper observations, 
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our hope is that it will be possible to discover evidence for sig- 
nificant gas infall in proto-galaxies, caught in the process of 
their initial assembly. 
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